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ABSTRACT 

I present the results of Monte-Carlo orbital simulations of Galactic Neutron Stars (NSs). The simulations 
take into account the up-to-date observed NS space and velocity distributions at birth, and account for their 
formation rate. I simulate two populations of NSs. Objects in the first population were born in the Galactic 
disk at a constant rate, in the past 12 Gyr. Those in the second population were formed simultaneously 12 Gyr 
ago in the Galactic bulge. I assume that the NSs born in the Galactic disk comprise 40% of the total NS 
population. Since the initial velocity distribution of NSs is not well known, I run two sets of simulations, 
each containing 3 x 10'' simulated NSs. One set utilizes a bimodal initial velocity distribution and the other 
a unimodal initial velocity distribution, both are advocated based on pulsars observations. In light of recent 
observational results, I discuss the effect of dynamical heating by Galactic structure on NS space and velocity 
distributions and show it can be neglected. I present catalogue of simulated NS space and velocity vectors in 
the current epoch, and catalogue of positions, distances and proper motions of simulated NSs, relative to the 
Sun. Assuming there are 10^ NSs in the Galaxy, I find that in the solar neighborhood the density of NSs is 
about 2 — 4 X 10^^ pc^^, and their scale height is about 0.3 — 0.6kpc (depending on the adopted initial velocity 
distribution). These catalogue can be used to test the hypothesis that some radio transients are related to these 
objects. 

Subject headings: stars: neutron — Galaxy: kinematics and dynamics 



1. INTRODUCTION 

The present day space and velocity distributions of Galac- 
tic isolated old Neutron Stars (NS), are the subject of many 
studies. Among the reasons for deriving the positional and 
kinematical properties of isolated old NSs was the suggestion 
that Gamma-Ray Bursts (GRBs) may originate from Galac- 
tic NSs (e.g., Mazets et al. 1980). This was refuted^ by the 
homogeneous sky distribution of GRBs (e.g., Meegan et al. 
1992), and later on by the discovery of their cosmic origin 
(e.g., Metzger et al. 1997; van Paradijs et al. 1997). 

Another exciting possibility, that re-ignited these efforts, 
was to detect isolated old NSs in soft X-ray radiation that may 
be emitted as they slowly accrete matter from the interstel- 
lar medium (ISM; Ostriker, Rees & Silk 1970; Shvartsman 
1971). Treves & Colpi (1991), and Blaes & Madau (1993) 
predicted that ^ 10^-10"*^ isolated old NS accreting from the 
ISM would be detected by the TJO^Ar Position Sensitive Pro- 
portional Counters (PSPC) all sky survey (Voges et al. 1999). 
Although intensive searches for these objects were carried out 
(e.g., Motch et al. 1997; Maoz, Ofek, & Shemi 1997; Haberl, 
Motch, & Pietsch 1998; Rutledge et al. 2003; Agueros et al. 
2006), only a handful of candidates were found. However, 
these are presumably young cooling NSs, rather than old NSs 
whose luminosity is dominated by accretion from the ISM 
(e.g., Neuhauser & Triimper 1999; Popov et al. 2000; Treves 
et al. 2001). Apparently, the reasons for the rareness of these 
objects in the ROSAT source catalog is that their typical ve- 
locities were underestimated by earlier studies (i.e., Narayan 
& Ostriker 1990). It is also possible that their accretion rate 
is below the Bondi-Hoyle rate^ (e.g., Colpi et al. 1998; Perna 
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- At least in the sense that the majority of GRBs are extragalactic; see 
Kasliwal et al. (2008). 
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et al. 2003; Toropina et al. 2001, 2003, 2005; see however 
Arons & Lea 1976,1980), or that these NSs are in an ejection 
stage (e.g., Colpi et al. 1998; Livio et al. 1998; Popov & 
Prokhorov 2000). 

During the last several years, new classes of radio transients 
were found (e.g.. Bower et al. 2007; Matsumura et al. 2007; 
Niinuma et al. 2007; Kida et al. 2008; see also Levinson et al. 
2002; Gal- Yam etal. 2006). Bower etal. (2007) found several 
examples of these transients in a single small field of view, 
and showed that they have time scales above 20 min but below 
seven days, and they lack any quiescent X-ray, visible-light, 
near-Infrared (IR), and radio counterparts. Recently, Ofek et 
al. (2009) suggested that these transient events may be asso- 
ciated with Galactic isolated old NSs. However, testing this 
hypothesis requires knowledge of the NSs space distribution. 

There are two approaches in the literature for calculating 
the theoretical space and velocity distributions of old NSs. 
Paczynski (1990), Hartmann, Woosley, & Epstein (1990), 
Blaes & Rajagopal (1991), Blaes & Madau (1993), and Pos- 
selt et al. (2008) carried out Monte-Carlo simulations of NS 
orbits. In such simulations, the positions and velocities of NS 
at birth are integrated assuming some non-evolving Galactic 
potential. A second approach is to use some sort of semi- 
analytic approximation in order to estimate the "final" NS 
space and velocity distributions. Frei, Huang, & Paczyn- 
ski (1992) calculated the final vertical density and velocity 
distributions of NSs, assuming that the gravitational poten- 
tial is a function only of the height, z, above the Galactic 
plane. Using the epicyclic approximation, Blaes & Rajagopal 
(1991) developed a technique that allows calculation of the 
full three dimensional velocity distribution of NSs. However, 
they showed that this method is not adequate for fast moving 
objects, which constitute the majority of the NS population. 
Blaes & Madau (1993) used the thin-disk approximation to 
calculate the space and velocity distributions of NSs. In this 
prescription, the radial motion of NSs are controlled by the 
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Galactic potential in the Galactic disk, regardless of the ver- 
tical height, z, of the NS above/below the Galactic plane. Fi- 
nally, another solution was presented by Prokhorov & Postnov 
(1994), who assume that the ergodic hypothesis is correct (see 
discussion in Binney & Tremaine 1987). 

Madau & Blaes (1994) noted that all these approaches ne- 
glect dynamical heating of NSs due to encounters with giant 
molecular clouds, spiral arms, and stellar "collisions" (e.g., 
Kamahori & Fujimoto 1986, 1987; Bai-banis & Woltjer 1967; 
Carlberg & Sellwood 1985; Jenkins & Binney 1990). They 
crudely estimated the order of magnitude of this effect by ap- 
plying the force-free diffusion equation to the vertical height, 
z, and NS speed distributions. They found that, at the so- 
lar neighborhood, dynamical heating may decrease the local 
density of NSs by a factor of ^ 1.5 relative to a non-heated 
population. 

In this paper 1 present the results of Monte-Carlo orbital 
simulations of Galactic NSs. These simulations improve 
upon past efforts (e.g., Paczynski 1990; Hartmann, Epstein, 
& Woosley 1990, Blaes & Rajagopal 1991; Blaes & Madau 
1993; Popov et al. 2005; and Posselt et al. 2008) in several as- 
pects. First, 1 use up-to-date space and velocity distributions 
of NSs at birth from Arzoumanian et al. (2002) and Faucher- 
Giguere & Kaspi (2006). Second, 1 assume a birth rate of 
NSs along the Galaxy life time, instead of assuming that all 
the NSs were born about 10 Gyr ago. Third, 1 generate a large 
sample of simulated NSs, which is about one to two orders of 
magnitude larger than those presented by previous efforts. 

The structure of this paper is as follows: In fJDl present the 
model ingredients for the Monte-Carlo simulations. In Sj3]l 
discuss dynamical heating and show that it can be neglected. 
The Monte-Carlo simulations are described in [jH and the cat- 
alogue of simulated NSs are presented in S|5] The results of 
these simulations are presented in !j6] and in i|2]l summarize 
the results and compare them with some of the previous ef- 
forts. 

2. MODEL INGREDIENTS 

In the following subsections, 1 present the ingredients of th e 
Monte-Carlo simulations. These are: NSs birt h rat e (i i2.1l i; 
space and velocity distributions of NSs at birth { ^2.2\ \ and the 
Galactic potential (i )2.3b . 

2.1. NSs birth rate 

The stellar population in the Milky Way is composed of 
at least two major components: a bulge and a disk. I there- 
fore, simulate two populations of NSs. A Galactic disk-born 
population and a Galactic bulge population. For the Galactic- 
disk NSs, 1 assume a continuous constant formation rate in the 
past 12 Gyr. This assumption is motivated by the analysis of 
Rocha-Pinto et al. (2000) that did not find any major trends in 
the star formation rate in the disk of our Galaxy. 1 assume that 
the Galactic -bulge NSs were born in a single burst 12 Gyr ago 
(Ferreras, Wyse, & Silk 2003; Ballero et al. 2007; Minniti & 
ZoccaH 2008). 

As I discuss in this section, the simulations presented here 
consist of 40% disk-born NSs and 60% bulge-born NSs. The 
actual ratio of disk to bulge NSs is unknown and may be sig- 
nificantly different than the one assumed here. Therefore, in 
the resulting catalogue (!|5]l I specify the origin of each simu- 
lated NS and present the results also for each population sep- 
arately. 

The total number of NSs in the Milky Way is constrained by 
the chemical composition of the Galaxy. Specifically, the iron 



content of the Galaxy implies that the total number of core 
collapse supernovae (SNe; and therefore NSs) that exploded 
in the Milky Way is about 10^ (Arnett, Schramm, & Truran 
1989). However, the star formation rate in the Galactic disk 
was approximately constant during the last ~ 12 Gyr (Rocha- 
Pinto et al. 2000; see also Noh & Scalo 1990). Therefore, 
with the current NS birth rate of up to one in 30 years (e.g., 
Diehl et al. 2006), one expects that there will be 4 x 10^ 
disk-born NSs in our Galaxy. The predictions based on the 
SN rate (< 4 x 10^) and the chemical evolution (^ 10^) are 
therefore inconsistent. Ways around this problem include un- 
certainty in the SN rate (e.g., Arnett et al. 1989), and difficulty 
in estimating the star formation history of our Galaxy (Keane 
& Kramer 2008). Another viable resolution of this discrep- 
ancy is that the Galactic bulge had a higher star formation 
rate at earlier times. 

The Galactic bulge contains a considerable number of stars 
that were born, apparently, in a ^ 1 Gyr-long, burst about 
12 Gyr ago (Ferreras, Wyse, & Silk 2003; Ballero et al. 2007; 
Minniti & Zoccali 2008). I therefore, assume here that in ad- 
dition to the disk-born NSs there are up to 6 x 10^ NSs that 
were formed in the Galactic bulge about 12 Gyr ago. This 
number was selected such that the total number of NSs in the 
Galaxy is lO**. 

1 note that the Galactic bulge contains only about 20% of 
the Galactic stars (e.g., Klypin et al. 2002). However, Ballero 
et al. (2007) argued that the initial mass function, dN /dM, of 
the bulge is skewed toward high masses, with x, the power- 
law index in the stellar mass function (dN /dM °^ M^''+'^^; 
e.g., Salpeter 1955), being < 0.95. Assuming that stars with 
masses exceeding 8 M,:: produce NSs within a relatively short 
time after their birth, this suggests that the number of NSs 
born in the Galactic bulge is a few times larger than the num- 
ber of NSs born in the Galactic disk, per unit stellar mass. 
Therefore, the assumption that the majority of NSs were born 
in the Galactic bulge is in rough agreement with the expected 
fraction of high mass (^ 8M0) stars in the bulge, relative to 
the disk. Nevertheless, as stated before, the ratio between the 
numbers of bulge-born to disk-born NSs is uncertain. The 
simulations presented here cover a wide range of possibilities 
regarding birth place and time scenarios. 

2.2. NSs space and velocity distributions at birth 

Following, Paczynski (1990), I adopt a radial birth proba- 
bility distribution, p{R), which follows the exponential stellar 
disk: 

dp{R)=aR-^e.x^[- — ]dR, (1) 

where R is the distance from the Galactic center projected on 
the Galactic plane, R^i — 4.5 kpc is the exponential disk scale 
length, and is the normalization given by: 

aR = [1 +R„^^^/R^)]-\ (2) 

where = 15 kpc, is the disk truncation radius. Yusi- 
fov & Ku9uk (2004) estimated the radial distribution (mea- 
sured from the Galactic center) of Galactic pulsars, and found 
it to be more concentrated than predicted by the Paczynski 
(1990) distribution. However, I note that the supernova rem- 
nant Galactic radial distribution is in good agreement with the 
Paczynski (1990) model. 

For the NS bulge population I assume that the space distri- 
bution at birth has the same functional form and parameters 
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as the NS disk population, but with R^i = 1 kpc in Eqs.[T]and 

El 

For the initial velocity probability distribution, I use two 
different models. Both models are consistent with the ob- 
served velocity distribution of young pulsars. One is a bi- 
modal velocity distribution composed of two Gaussians found 
by Arzoumanian et al. (2002), and the second is based 
on a double-sided exponential unimodal velocity distribution 
(Faucher-Giguere & Kaspi 2006). 

Bimodal velocity distribution: The first initial velocity 
distribution we used is from Arzoumanian et al. (2002). It 
consists of two components. The three-dimensional speed, v, 
probability density of this model is: 

Mv) ^ W^g^ (24i)^/^ exp(-vV[2cT,?,]) (3) 

where j G {1,2}, /i = 0.40 and /2 = 0.60 are the fractions 
of the two NS velocity components, C7, i = 90kms^', and 
(7,,2 =500kms-i. 

Unimodal velocity distribution: The second model I use 
for the velocity distribution of NSs at birth was advocated by 
Faucher-Giguere & Kaspi (2006). In their model, each one- 
dimensional (projected) velocity component, v,, consists of a 
double-sided exponential probability distribution: 

p(v,) = -^exp(-M), (4) 

2{vi) \ {vi)J 

where (v,) = 180km s^' is the mean velocity in each direc- 
tion, /. I note that drawing the velocity components of simu- 
lated NSs directly from Eq.|4]will result in an aspherical dis- 
tribution. In order to avoid this problem the corresponding 
three-dimensional distribution is required. Faucher-Giguere 
& Kaspi (2006) noted that this distribution is difficult to derive 
analytically, and did not provide its functional form. How- 
ever, it is possible to show that the corresponding spheri- 
cally symmetric three-dimensional distribution corresponding 
to Eq. |4]is given by 

^(,) = _J_exp(-|i). (5) 

In Figure [T] I show the one-dimensional (dashed lines) and 
three-dimensional {solid lines) initial velocity probability dis- 
tributions based on Faucher-Giguere & Kaspi (2006) veloc- 
ity distribution (thin black lines) and the Arzoumanian et al. 
(2002) distribution {thick gray lines). Several other fits for 
the pulsars velocity distribution at birth are available (e.g., 
Hansen & Phinney 1997; Cordes & Chernoff (1998); Hobbs 
et al. 2005; Zou et al. 2005). However, the Arzoumanian 
et al. (2002) and Faucher-Giguere & Kaspi (2006) velocity 
distributions cover a wide range of possibilities, that reflect 
the current uncertainties in the pulsars velocity distribution at 
birth. 

The two models are considerably different from each other, 
mainly in the low and high velocity tails. Relative to the 
Arzoumanian et al. (2002) velocity distribution the Faucher- 
Giguere & Kaspi (2006) distribution has excess of NSs at low 
velocities. Since the slowest NSs are the most efficient accre- 
tors from the ISM (e.g., Ostriker et al. 1970), these differences 
may have a significant impact on the detectability of isolated 
old NSs in X-ray wavebands (e.g., Blaes & Madau 1993). 

The velocity distributions, in Eqs.|3]and|5] are given relative 
to a rotating disk. Therefore, I added the disk rotation velocity 
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Fig. 1. — The initial velocity distribution of NSs at birth relative to a ro- 
tating galaxy. The dashed lines represent the one-dimensional (ID) velocity 
distributions, while the solid lines show the three-dimensional (3D) speed dis- 
tributions. The black thin lines are for the Faucher-Giguere & Kaspi (2006) 
unimodal double-sided exponential velocity distribution, while the gray thick 
lines show the Arzoumanian et al. (2002) bimodal Gaussian velocity distri- 
bution. 

to these birth velocities 

v',^V,-Vcirc{R)MQ). (6) 

v;.-vv + y„>,(/?)cos(0), (7) 

v; = v„ (8) 

where v, is the velocity vector relative to the rotating disk (i.e., 
the one obtained from Eqs.[3]or|5l), vj is the velocity vector in 
a non-rotating (inertial) reference frame, R = (;ic^ 6 is 

azimuthal angle on the Galactic disk {— arctan2[3',z]), which 
was selected from a uniform random distribution, and x, y, 
and z are the position in a Cartesian, non-rotating, coordinate 
system, whose origin is the Galactic center, and z is vertical 
to the Galactic plane. Finally, Vdrc was obtained from 

Vdrc{R) = V^R^, (9) 

where <I> is the Galactic potential (see il2.3l l. at the point of 
interest. 

Finally, I assume that the vertical height above the Galactic 
mid-plane of NSs at birth, z, is drawn from a Gaussian distri- 
bution (o= exp[— z^/ (2aft)]), with at — 0.16 kpc and 0.05 kpc, 
for the Arzoumanian et al. (2002) and Faucher-Giguere & 
Kaspi (2006) initial velocity distributions, respectively. 

2.3. The Galactic potential 

Following Paczynski (1990), I assume that the Galactic po- 
tential is composed of a disk component, an spheroidal com- 
ponent, and an halo component. I also added a central black 
hole component. 

The disk and spheroid components are described by the fol- 
lowing potential proposed by Miyamoto & Nagai (1975), 

^kix,y,z) = — ^ -jj^, (10) 

[x^+y^ + [ak+{z^ + bly/2]^j 

where a/., and M/. are given in Equations [T3] and [141 and 
where k = d corresponds to the disk component and k — s 
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corresponds to the spheroid component, and G is the Gravita- 
tional constant. Next, the halo potential is given by 



GM,, 



were — y? 



,.2 I 



1 



■Inl 



■ arctan 



(11) 



Further- 



y"+z^, Mij, and r/, are listed in Eq.| 
more, I added a component representing the Galactic central 
massive black-hole (e.g., Ghez et al. 1998): 



GMhh 



(12) 



where M/^/, is given in Eq. [16] (Eisenhauer et al. 2005). 
The choice of parameters listed here reproduces the observed 
Galactic rotation, local density, and local column density (see 
Paczynski 1990 for details): 



a,i = 3.7 kpc, b,i = 0.20 kpc, M,, = 8.07 x 10^" M^iU) 



= O.Okpc, /7, = 0.277 kpc, M, = 1.12 x 10^° M0(14) 




Fig. 2. — Stellar velocity dispersion (a) along the U (radial; circles), V 
(tangential; triangles) and W (vertical; squares) Galactic directions as a func- 
tion of age since birth, T, as mea sured by Nordstrom et al. (2004). The lines 
represent the best fit of Equation [T8]to the data. The best fit parameters are 
listed in Table[T] 



r,, = 6.0 kpc, M/, = 5.0 x 10^" Mq, (15) 

Mbh = 3.6 X 10^ Mq. (16) 

I note that the black hole contributes about 4% (0.6%) of 
the gravitational potential at distance of 1 pc (lOpc) from the 
black hole. Therefore, its influence on the Galactic potential 
is negligible and it does not change the fitted parameters in 
Eqs. [T3l[T6] However, it may heat NSs passing nearby. Fi- 
nally, the Galactic potential is the sum of these four compo- 
nents: 

* = 4>,, + 4>, + 4>A + *w,. (17) 

3. DYNAMICAL HEATING 

Dynamical heating, presumably by giant molecular clouds 
(e.g., Kamahori & Fujimoto 1986; 1987), spiral structure 
(e.g., Barbanis & Woltjer 1967; Carlberg & Sellwood 1985; 
Jenkins & Binney 1990), and stars, tends to broaden the ve- 
locity and spatial distributions of Galactic stars (e.g., Wielen 
1977; Nordstrom et al. 2004). 

In order to roughly estimate the effect of dynamical heat- 
ing on NSs, Madau & Blaes (1994) applied the force-free 
diffusion equation to the velocity and vertical distance dis- 
tributions of NSs. They adopted the total diffusion coef- 
ficient, C (= Cu +Cv +Cw), measured by Wielen (1977; 
C = 600km^ s^^ Gyr^'). In this, U corresponds to the ra- 
dial direction, and is positive in the direction of the Galactic 
center; V points in the direction of circular rotation; and W is 
directed towards the North Galactic pole. Applying heating, 
they found that the local density of NSs is smaller by about 
30% relative to the case of no heating. In the following we 
discuss this approach in light of new observational data avail- 
able. 

We estimate the diffusion coefficient, C, using modern data. 
Nordstrom et al. (2004), estimated the dynamical heating by 
the Galactic disk and showed that it does not saturate after 
some time as suggested by measurements based on smaller 
samples (e.g., Quillen & Garnett 2001; see also Aumer& Bin- 
ney 2009). Approximating dynamical heating by a random 



TABLE 1 

The Galactic disk diffusion coefficients 



Velocity component 




C 


X'/dof 




km s ' 


km^ s-2 Gyr-' 




U 


18.7±1.2 


185 ±19 


3.1/8 


V 


9.7 ±0.9 


81.7±7.5 


9.3/8 


w 


5.5±1.0 


61.3±5.1 


11.0/8 



Note. — The diffusion coefficients of Galactic stai's as obtained by fitting 
Eq. llSl with the measurements of Nordstrom et al. (2004). 



walk process, I fit the Nordstrom et al. (2004) measurements 
with the function (e.g., Wielen 1977) 

a,^{alo + C,t)'/\ (18) 

where CT^, is the velocity dispersion component at time, t, 
since birth, a^ o is the initial velocity dispersion component 
and g G {U,V,W). The data and the best fit curves are shown 
in Figure |2] and the best fit parameters are listed in Table [T] 
The total diffusion coefficient that I find, 
Ct/ + Cv + Cw = 328 ± 2 1 km^ s^^ Qyy-- 1 ^ about half of the 
value found by Wielen (1977), and used by Madau & Blaes 
(1994). 

Next, the initial velocities estimated by Narayan & Ostriker 
(1990), that was used by Madau & Blaes (1994), are consider- 
ably lower than the more recent estimates (e.g., Arzoumanian 
et al. 2002). As suggested by Eq.[T8] diffusion affects mostly 
low velocity objects. Therefore, dynamical heating is less im- 
portant than estimated by Madau & Blaes (1994). 

Furthermore, the approach taken by Madau & Blaes (1994) 
assumes that NSs are being affected by diffusion at all times. 
However, the scatterers are restricted to the Galactic plane. 
NSs are born with high velocities and spend 80% to 90% 
of the time at distances larger than lOOpc from the Galactic 
plane (based on the results in S]6]l. Therefore, they are less sus- 
ceptible to dynamical heating than disk stars. Based on these 
three arguments, I conclude that the importance of dynamical 
heating was over estimated by Madau & Blaes (1994) by at 
least an order of a magnitude. 
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Nevertheless, dynamical heating may effect some of the 
slow moving objects, but they are a minority among the 
Galactic NSs. Eq.fTSlroughlv suggests that dynamical heating 
is important for NSs with speeds smaller than about (Ct)'/^. 
Even for T = lOGyr this gives 60kms^^ However, only 
2.8% and 4.5% of the NSs in the Arzoumanian et al. (2002) 
and Faucher-Giguere & Kaspi (2006), initial velocity distri- 
butions, respectively, have speeds smaller than 60kms^^ at 
birth (relative to their local standard of rest). 

4. MONTE-CARLO SIMULATIONS 
To solve for NS orbits I integrate the equations of motion 

(fix ^ 

^ = -V*, (19) 

using a Livermore ordinary differential equation solver^ 
(Hindmarsh 1983). The integration is performed in a non- 
rotating, Cartesian coordinates system, the origin of which 
is the Galactic center, and the z axis is perpendicular to the 
Galactic plane. 

For each initial velocity distribution (i.e., bimodal or uni- 
modal; see il2.2l ). I simulated 3 x lO*" NS orbits. In each simu- 
lation I randomly drew the NSs birth times, positions, an d ve- 
lociti es, f rom the probability distributions described in i ]2.1l 
and i l2.2l As explained before, 40% of these NSs are disk- 
born, and the rest are bulge-born. 

At the end of each simulation I checked if the integration 
conserved the total energy. In cases in which the energy was 
not conserved to within 0. 1 %, I reran the integration using the 
same initial conditions, with a refined integration tolerance. In 
the second (and final) iteration, the energy was conserved to 
better than 2.8% in all cases. 

5. CATALOGUE 

The catalogue of initial (i.e., at birth) and final (i.e., cur- 
rent epoch) simulated NSs space and velocity components are 
listed in Table |2] for the bimodal initial velocity distribution 
of Arzoumanian et al. (2002), and in Table [3] for the uni- 
modal initial velocity of Faucher-Giguere & Kaspi (2006). 
The first column in each table indicates whether the simulated 
NS belongs to the bulge population (code 0) or disk popula- 
tion (code 1). 

In addition, catalogue of simulated NS positions, distances, 
radial velocities, and proper motions for an observer located 
at the "solar circle", moving around the Galactic center with a 
velocity of 220 km s^^ is given in tables|4]and|5] The solar cir- 
cle is defined to be on the Galactic plane (i.e., z = kpc) at the 
distance of the Sun from the Galactic center {Re, = 8.0 kpc; 
Ghez et al. 2008). The catalogue in Tables liandlS] are 
based on the initial velocity distribution of Arzoumanian et 
al. (2002) and the Faucher-Giguere & Kaspi (2006), respec- 
tively. I note that the radial velocity and proper motions are 
calculated for a static observer with respect to the Local Stan- 
dard of Rest (LSR; i.e., the solar motion with respect to the 
LSR is neglected). 

In order to reduce Poisson errors in the local properties of 
NSs (i.e., density; sky surface density), tables |4] and |5] were 
produced by calculating the positions of the 3x10^ simulated 
NSs in tables |2]and|3] respectively, from 100 random locations 
on the solar circle. Hence, tables|4]and|5]list 3x10^ simulated 
NSs. 

^ http://www.netlib.org/odepack/ 
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Fig. 3. — The total energy distribution of Galactic-bom NSs. The energy 
does not include external potentials (i.e., other galaxies), and assumes all NSs 
have a mass of 1 .4 M j, . The different lines shows the distribution of all the 
NSs (solid thick Mack lines); only disk-born NSs {solid thin lines); and only 
bulge-bom NSs (dashed lines). 

6. STATISTICAL PROPERTIES 

In this section I present the space and velocity distributions 
of simulated NSs, at the current epoch. In il6.1I I discuss the 
overall distribution of NSs in the Galaxy, while in il6.2| I dis- 
cuss their statistical properties as observed from the LSR. Ad- 
ditional specific statistical properties of these objects are dis- 
cussed in Ofek et al. (2009) in the context of the long-duration 
radio transients (Bower et al. 2007; Kida et al. 2008). The re- 
sults presented in this section assume there are N<) = N 
NSs in the Galaxy, where is the number of NSs in the 
Galaxy, of which 40% were born in the disk and 60% in the 
bulge. 

6.1. Overall properties 

NSs are born with large space velocities, which are typ- 
ically of the order of the escape velocity from the Galaxy. 
These are presumably the result of kick velocities due to 
asymmetric supernovae explosions (e.g., Blaauw 1961 ; Lai 
et al. 2006). Therefore, it is expected that a large fraction of 
the Galactic NSs will be unbounded to the Milky Way grav- 
itational potential and some may be found at very large dis- 
tances. 

In Figure |3] I show the distribution of the total (kinetic 
+ potential) energy of the simulated NSs: Ma?5(v^/2 + <I>), 
where the NSs mass was set to M[^s = 1-4M0. Panel (a) 
shows the energy distribution based on the simulations using 
the unimodal initial velocity distribution of Faucher-Giguere 
& Kaspi (2006), while panel (b) is for the bimodal distribu- 
tion of Arzoumanian et al. (2002). In each panel, the thick 
black solid line represents the entire NS population, the thin 
solid line shows the disk-born NSs, and the dashed line is 
for bulge-born NSs. Using the approximation that all the 
NSs with negative energies are gravitationally bounded to the 
Galaxy (neglecting heating), in Table|6]l give the fractions of 
NSs bounded to the Galactic gravitational potential. 

Given the large fraction of gravitationally-unbounded NSs, 
it is expected that some NSs may be found at very large dis- 
tances, r, from the Galactic center. I find that about 12% 
and 35% of the NSs born in the Galaxy are currently at dis- 
tances larger than 1 Mpc from the Galactic center, for the 
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TABLE 2 

NS Monte-Carlo orbital simulations using the bimodal initial velocity distribution 













Initial 














Final 






p 


Age 


X 


y 


Z 


X 


Y 


Z 




X 


Y 


Z 


X 


y 


Z 




Gyr 


kpc 


kpc 


kpc 


kpc Gyr^ ' 


kpc Gyr" ' 


kpc Gyr" ' 


kpc Gyr" ' 


kpc 


kpc 


kpc 


kpc Gyr" ' 


kpc Gyr"' 


kpc Gyr" ' 





12.0 


-2.04 


1.16 


-0.07 


-11.1 


-194.6 


8.8 


208.7 


-0.71 


2.41 


-0.07 


-181.0 


-70.7 


9.1 





12.0 


-0.58 


1.57 


-0.27 


-199.1 


0.1 


-376.7 


206.7 


-1.40 


-8.23 


-5.08 


42.4 


25.4 


-70.3 





12.0 


0.60 


1.57 


0.02 


531.9 


507.1 


-60.9 


206.6 


477.31 


404.37 


-5.35 


-154.9 


-132.3 


2.3 





12.0 


-1.27 


2.33 


0.05 


-140.6 


-96.1 


-53.5 


211.6 


-0.20 


2.64 


0.06 


-167.7 


-32.5 


52.4 





12.0 


-1.54 


-2.95 


-0.18 


-137.7 


-209.2 


-425.3 


218.2 


-1.31 


-1.29 


-19.70 


-10.4 


53.7 


251.2 



Note. — Catalog of initial and final space positions and velocity components for 3 X 10* simulated NSs. The NSs were simulated using the bimodal initial 
velocity distribution of Arzoumanian et al. (2002). The velocity components are given in kpc Gyr"', relative to the non-rotating Galaxy. In order to convert 
kpc Gyr"' to kms"', divide it by 1.0227. P is the population type: for bulge-born NS; 1 for disk-bom NS. Vdrc refers to the Galaxy rotation speed at the 
projected location, on the Galactic disk, in which the NS was born. The initial (and final) velocity components are given relative to a non-rotating galaxy (ineilial 
reference frame). The numbers in this table are rounded in order to fit into the page. This table is published in its entirety in the electronic edition of this paper. 
A portion of the full table is shown here for guidance regarding its form and content. 



TABLE 3 

NS Monte-Carlo orbital simulations using the unimodal initial velocity distribution 













Initial 














Final 






p 


Age 


X 


y 


Z 


X 


Y 


Z 




X 


y 


Z 


X 


Y 


Z 




Gyr 


kpc 


kpc 


kpc 


kpc Gyr" ' 


kpc Gyr" ' 


kpc Gyr" ' 


kpc Gyr" ' 


kpc 


kpc 


kpc 


kpc Gyr" ' 


kpc Gyr" ' 


kpc Gyr" ' 


1 


2.3 


-7.05 


-6.69 


0.07 


-258.0 


-232.1 


-55.2 


220.2 


-39.41 


-3.21 


-5.51 


24.9 


4.3 


8.8 





12.0 


0.38 


0.88 


-0.01 


-124.3 


-176.5 


134.8 


226.7 


0.35 


0.10 


0.60 


213.6 


178.4 


79.1 


1 


2.9 


-0.69 


5.40 


0.07 


-124.5 


122.9 


294.4 


227.1 


-6.97 


-5.47 


-11.00 


34.2 


-57.3 


105.9 


1 


0.6 


-2.02 


-1.49 


-0.02 


-510.5 


-374.6 


5.6 


210.2 


-140.81 


-103.14 


-0.48 


-139.0 


-101.7 


-0.8 





12.0 


-2.34 


-2.96 


-0.02 


178.0 


-374.0 


-108.7 


221.6 


-7.62 


-17.92 


2.54 


115.5 


87.9 


-21.8 



Note. — Like Table|2] but for the unimodal initial velocity distribution of Faucher-Giguere & Kaspi (2006). 



TABLE 4 

Catalog of simulated NSs as observed from the LSR - based on the bimodal initial velocity distribution 



p 


Age 


/ 


b 


dist 






RV 




Gyr 


deg 


deg 


kpc 


"yr"' 


"yr"' 


kms"' 





12.0000 


354.72497 


-0.39865 


10.0003 


-0.0064166 


0.0001654 


-149.7 





12.0000 


72.14501 


-38.10227 


8.2258 


-0.0025465 


-0.0041258 


-93.8 





12.0000 


273.07904 


-0.49003 


625.9548 


-0.0000573 


0.0000017 


340.7 





12.0000 


351.53160 


0.34229 


9.9375 


-0.0057996 


0.0011040 


-124.8 





12.0000 


13.81616 


-68.84354 


21.1283 


-0.0015999 


0.0004210 


-247.1 



Note. — Catalog of Galactic longitudes (/), latitudes (b), distances, proper motions in Galactic longitude (/i;) and latitude (^l;,), and radial velocities (RV) 
from a point on the solar circle (e.g., the LSR) of 3 x 10* simulated NSs. The velocities and proper motions do not include the motion of the Sun relative to the 
LSR. The proper motions are given in the Galactic coordinate system. The catalog was generated by calculating the positions of the 3 x 10* NSs in Table|2](i.e., 
assuming the Arzoumanian et al. [2002] initial velocity distribution) as observed from 100 random points on the solar circle. The conversion of space velocity 
to proper motion and radial velocity was carried out using the inverse of Eq. 3.23-3 in Seidelmann (1992, p. 121). This table is published in its entirety in the 
electronic edition of this paper. A portion of the full table is shown here for guidance regarding its form and content. 



TABLE 5 

Catalog of simulated NSs as observed from the LSR - based on the unimodal initial velocity distribution 



p 


Age 


/ 


b 


dist 




Mi, 


RV 




Gyr 


deg 


deg 


kpc 


"yr"' 


"yr"' 


kms"' 


1 


2.2753 


237.98450 


-8.95638 


35.3986 


0.0008045 


0.0002084 


166.7 





12.0000 


2.28569 


4.42334 


7.8361 


-0.0014499 


0.0016472 


212.2 


1 


2.8584 


306.10398 


-45.37962 


15.4574 


-0.0018509 


0.0033513 


96.8 


1 


0.6194 


262.62769 


-0.15823 


173.4066 


-0.0001141 


0.0000002 


334.3 





12.0000 


308.68736 


6.21614 


23.4368 


0.0000392 


-0.0003617 


171.9 



Note. — 



Like Table|5] but for the unimodal initial velocity distribution of Faucher-Giguere & Kaspi (2006). 
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TABLE 6 

Fraction of NSs gravitationally bounded to the Galaxy 



Initial velocity" 


All 


disk 


bulge 


A2002 


0.38 


0.16 


0.52 


FK2006 


0.30 


0.13 


0.41 



" Initial velocity distribution used in the simulations, where A2002 corre- 
sponds to Arzoumanian et al. (2002), and FK2006 to Faucher-Giguere & 
Kaspi (2006). 

unimodal and bimodal initial velocity distributions, respec- 
tively. For 1 Mpc< r < lOMpc, I find that the density of NSs 



1.9 X 10 

-3 



-5 



(r/lkpc)-^'^A/9pc 



as a function of r is about: 

and 3.6 x 10^^(r/lkpc)^^ 'A'9 pc^^, for the bimodal and uni- 
modal initial velocity distributions, respectively. Finally, I 
find that some Milky Way born NSs may be at distances as 
large as 30 to 40Mpc from the Galaxy. I note that the lo- 
cal density, in our Galaxy, of NSs born in other galaxies, is 
of the order of lO^^A'gpc^^. This was estimated by calcu- 
lating the density in the the Milky Way, of NSs born in each 
galaxy found within lOMpc. For this, I used a version of the 
Tully (1988) nearby galaxy catalog (Ofek 2007) where the to- 
tal number of NSs in each galaxy was normalized by its total 
Z?-band magnitude, relative to Milky Way. 

In Figures m and |5] I show the density of NSs on the Galac- 
tic plane, and the surface density of NSs projected on the 
Galactic plane, respectively, as a function of distance from 
the Galactic center The notations are the same as in Figure[3] 
In addition, in Fig. |5] the gray solid line represents the ini- 
tial surface distribution of all NSs. Furthermore, the vertical 
dashed Unes mark the distance of the Sun from the Galactic 
center, Rq = S.Okpc (Ghez et al. 2008). As noted before, I 
refer to this distance from the Galactic center, when located 
on the Galactic plane, as the solar circle. The densities in the 
solar circle and Galactic center are listed in the figure cap- 
tions. I note that the fluctuations seen in small and large radii 
are due to Poisson noise. 

In Figure |6] I show the NSs density as a function of height 
above or below the Galactic plane, as measured at the solar 
circle. This was calculated by counting the number of simu- 
lated NSs with distance from the Galactic center, projected on 
the Galactic plane, of 8.0±0.5kpc. The line scheme is the 
same as in Figures [3] At the solar circle the scale height^ of 
NSs is about 0.6 kpc and 0.3 kpc for the Arzoumanian et al. 
(2002) and Faucher-Giguere & Kaspi (2006) initial velocity 
distributions, respectively. 

Shown in Fig.|2]are the space densities of NSs as a function 
of projected (on the Galactic plane) distance from the Galac- 
tic center, for several different Galactic height: kpc; 1 kpc; 
3 kpc; and 10 kpc. The densities are calculated in slices, par- 
allel to the Galactic plane, with semi-width of 0. 1 kpc. 

Finally, in Figure [8] I show the initial and final speed dis- 
tributions as measured relative to an inertial reference frame 
(contrary to Fig. [T] which shows the initial speed distribution 
relative to a rotating reference frame). The probabilities in 
this Figure are shown per Ikms^' bins. The line scheme 
is again like the one used in Fig.|5] As expected, the typical 
speeds of NSs decrease with time as they, on average, increase 
their distances from the Galaxy and lose kinetic energy. For 
the unimodal initial velocity distribution, at the current epoch 
I estimate that about 17%, 53%, and 99% of the NSs have 
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Fig. 4. — Density of NSs at the Galactic plane, as measured at the cun'ent 
epoch. See Fig.|3]for details regarding line types. The density is calculated 
by counting the NSs within 100 pc from the Galactic plane, dividing it by 
the appropriate volume in each bin, and assuming that the total number of 
NSs in the Galaxy is lO'. For the initial velocity distribution of Faucher- 
Giguere & Kaspi (2006; panel a), the NSs density at the solar circle (distance 
of S.Okpc from the Galactic center; Ghez et al. 2008) is 4.0 x lO^^Wg pc"^ 
of which 60% are disk born and 40% are bulge bom. The NSs density at the 
Galactic center is: 3 x 10^' Ng pc^^, of which 5% are disk born and 95% are 
bulge bom. For the initial velocity distribution of Arzoumanian et al. (2002; 
panel b), the density at the solar circle is 2.4 x lO^'^A'ij pc^^, of which 66% 
are disk born and 34% are bulge bom. The density at the Galactic center is: 
2 X lO^^Ng pc^^, of which 7% are disk born and 93% are bulge born. I note 
that the fluctuations seen in small and large radii are due to Poisson noise. 



10' 



,10" 



10" 



(a: Unimodal) 




(b: Bimodal) 




10 10 10 10 

R [kpc] 



10 10 10 10 10 

R [kpc] 



Fig. 5. — The surface density of NSs projected on the Galactic plane, as 
measured at the current epoch. See Fig.|3]for details. In addition the gray 
lines show the initial suii'ace density distributions. For the initial velocity 
distribution of Faucher-Giguere & Kaspi (2006; panel a), the NSs surface 
density at the solar circle is 0.6Nq pc^^ of which 54% are disk born NSs, 



The surface density at the Galactic center is 
5 are disk bom NSs, and the rest are bulge- 



and 46% are bulge-bom NSs. 
~ 90^9 pc^^ of which about ' 
born NSs. For the initial velocity distribution of Arzoumanian et al. (2002; 
panel b), the NSs density at the solar circle is OANg pc"- of which 61% are 
disk born NSs, and 39% are bulge-born NSs. The coiresponding surface 
density at the Galactic center is ~ 50^9 pc"- of which ~ 6% are disk bom 
NSs, and the rest are bulge-born NSs. As before, the plotted densities assume 
the total number of NSs is lO'. 



' Scale height is defined as the height at which the density drops by l/e. 
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Fig. 6. — The density of NSs as a function of heiglit, 7„ above/below tlie 
Galactic plane as measured on the solar circle at the current epoch. Line 
types are like those in Fig. [5] For the initial velocity distribution of Faucher- 
Giguere & Kaspi (2006; panel a), 50% (90%) of the NSs are found within 
0.9 kpc (5.6 kpc) from the Galactic plane. For the initial velocity distribution 
of Arzoumanian et al. (2002; panel b), 50% (90%) of the NSs are found 
within 0.7 kpc (4.4 kpc) from the Galactic plane. 
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Fig. 8. — The initial and final speed probability distributions of the simu- 
lated NSs, relative to an inertial reference frame. Line types are like those in 
Fig. [3] The probabilities are calculated for 1 kms"' bins. Panel (a) shows the 
results for the unimodal initial velocity distribution, while panel (b) for the 
bimodal initial velocity distribution. 
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Fig. 7. — The density of NSs in Galactic heights of Okpc, 1 kpc, 3 kpc, 
and 10 kpc above/below the Galactic plane, as a function of the projected (on 
the Galactic plane) distance from the Galactic center, R, and as measured at 
the current epoch. Panel (a) shows the densities calculated using the initial 
velocity distribution of Faucher-Giguere & Kaspi (2006). For this model, at 
the solar circle, the densities are about 4 X IQ-'^Nt) pc"^ 6 x IQ-^Nt) pc"^ 

2 X 10"^A/9 pc"^ and 2 x W^Nq pc"^, at Galactic heights of Okpc, 1 kpc, 

3 kpc, and 10 kpc, respectively. Panel (b) shows the densities calculated using 
the initial velocity distribution of Arzoumanian et al. (2002). For this model, 
at the solar circle, the densities are about 2 X 10"" pc"^ 5 x lO^^/Vg pc"-\ 
1 X lO^^A'gpc"^, and 1 X lO^^Wg pc"', at Galactic heights of Okpc, Ikpc, 
3 kpc, and 10 kpc, respectively. 

speeds below 100, 200, and lOOOkms^', respectively. For 
the bimodal initial velocity distribution, at the current epoch 
I estimate that about 11%, 39%, and 94% of the NSs have 
speeds below 100, 200, and 1000km s"\ respectively. 

6.2. Properties observed from the LSR 

In this subsection I discuss: (i) the expectancy distance of 
the nearest NS to the Sun; (ii) the number of young NSs in 
the solar neighborhood; (iii) the proper motion distribution of 
nearby NSs; and (iv) the all-sky distribution of Galactic NSs. 



The probability to find a NS within distance d from the Sun 
is given by 

P<rf-l-exp(-^7rpfl'^). (20) 

Given the local density of NSs, p, that I have found in ^6.\\ 
this implies that the expectancy distance of the nearest NS is 

— 1/3 —1/3 

about 8.8A^g pc and 7.5A^g pc for the Arzoumanian et 
al. (2002) and Faucher-Giguere & Kaspi (2006) distributions, 
respectively. 

For the bimodal initial velocity distribution (Table |4|, 
within Ikpc from the Sun, 63% of the NSs are disk-born, 
and there are about 220A^9 (900A^9) NSs younger than 1 Myr 
(lOMyr). On the other hand, for the unimodal initial velocity 
distribution (Table [S]), within 1 kpc from the Sun, 57% of the 
NSs are disk-born, and there ai-e about 190A^9 (930A^9) NSs 
younger than 1 Myr (10 Myr). 

In Figure |9l I show the median total proper motion {solid 
lines) of simulated NSs as a function of their distance from 
an observer located on the solar circle. The dotted lines show 
the lower and upper 95 -percentiles of the proper motion dis- 
tributions. The black lines represent the unimodal initial ve- 
locity distribution and the gray lines are for the bimodal ini- 
tial velocity distribution. In addition, the dots show the ob- 
served proper motions of known pulsars^ (Manchester et al. 
2005). I note that the largest total proper motion in Ta- 
ble H] (i.e., bimodal initial velocity distribution) is 4.6" yr"^ 
and that ^ nONg NSs are expected to have proper motion in 
excess of l"yr^'. In Table |5] the largest proper motion is 
3.3" yr"', and about 240A^9 NSs are expected to have proper 
motion in excess of 1" yr^'. 

In Fig. [Tok-f I show the sky surface density of NSs at the 
current epoch, for an observer located at the solar circle. Pan- 
els (a)-(c) are for the unimodal initial velocity distribution, 
while panels (d)-(f) are for the bimodal initial velocity distri- 
bution. Panels (a) and (d) show the distribution for all the NSs 

* Pulsar distances and proper motions obtained from: 
http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 9. — The median total proper motion (solid lines) of simulated NSs as 
a function of their distance from an observer located on the solar circle. The 
dotted lines show the lower and upper 95-percentiles of the proper motion 
distributions. The black lines represent the unimodal initial velocity distri- 
bution and the gray lines are for the bimodal initial velocity distribution. In 
addition, the dots show the observed proper motion of known pulsars. 

(i.e., disk- and bulge-born populations), panels (b) and (e) for 
the disk-born population, and panels (c) and (f) for the bulge- 
born population. The maps are shown in the Aitoff equal-area 
projection and Galactic coordinate system, where the Galac- 
tic center is at the center of each map. The surface densities 
are normalized assuming that there are 4x10^ disk-born NSs, 
and 6 x 10** bulge-born NSs. 

At the positions with the lowest surface density of NSs, the 
Poisson eiTors due to the limited statistics are smaller than 
about 10%. I find that the minimum surface density is at- 
tained at the direction of the Galactic poles, and is about 
3900A^9deg-2 and SlOOA^gdeg^^ for the unimodal and bi- 
modal initial velocity distributions, respectively. The maxi- 
mum surface density is at the direction of the Galactic cen- 
ter and it is about 1.4 x IQ^N<) deg"^ and 1.1 x lO^A^g deg"^, 
for the unimodal and bimodal initial velocity distributions, re- 
spectively. I note that the differences between the sky surface 
densities resulting from the two initial velocity distributions, 
are as large as about 65%. 

7. SUMMARY 

The Milky Way's NSs content, and in particular the NS 
space and velocity distributions are important for searching 
nearby isolated old NSs, and studying any ongoing activity 
from such objects. As a tool for such studies, I present a mock 
catalogue of simulated isolated old NSs spatial positions and 
velocities, at the current epoch. 

The catalogue were constructed by integrating the equa- 
tions of motion of simulated NSs in the Galactic potential. 
The simulations include two populations of NSs, one in which 
the NSs were born in the Galactic bulge about 12Gyr ago, 
and the second population in which NSs are being born at the 
Galactic disk at a constant rate, starting 12 Gyr ago. The com- 
bined NS population assumes that 60% of the NSs originated 



in the bulge and the rest in the disk. Although we do not know 
what is the exact number of Galactic NSs, and what was their 
position-dependent birth rate, these two populations provide 
a wide range of initial conditions. 
I generated two catalogue of simulated NSs. Each cata- 

TABLE 7 
Comparison with previous works 



pc^^ 




Z|/2 

kpc 


mode(/)[i'])" 
kms"' 


Reference 


2.4 X 10" 


4 


0.42 


240 


This Work (disk+bulge) A2002'' 


4.0 X 10" 


4 


0.20 


200 


This Work (disk+bulge) FK2006'= 


1.4 X 10" 


3 


0.20 




Paczynski (1990) 


7.5 X 10" 


4 


0.27 


43 


Blaes & Madau (1993) 


5.3 X 10" 


4 


0.50 


69 


Madau &Blaes (1994) 


4 X 10-"* 






140 


Pema et al. (2003) 


5 X 10"'' 




0.4 




Popov et al. (2005) 



Note. — The density, p, at the solar neighborhood is calculated assuming 
there are lO' NSs in the Galaxy. Zi^2 is the height above the Galactic plane, 
in which the NS density drops to 1/2 its value on the Galactic plane, calcu- 
lated at the solar circle. The density from Madau & Blaes (1994) is calculated 
assuming that diffusion operats for 5 Gyr (see [|3). The initial velocity distri- 
bution used by Paczynski (1990), Blaes & Madau (1993), and Madau & Blaes 
(1994) is from Narayan & Ostriker (1990). The initial velocity distribution 
in Pema et al. (2003) is from Cordes & Chemoff (1998). 
" Mode is the most probable value of the distribution, and the speeds are 
measured relative to an inertial reference frame. 

^ Using the bimodal initial velocity distribution of Ar'zoumanian et al. (2002). 

Using the unimodal double-sided exponential initial velocity distribution of 
Faucher-Giguere & Kaspi (2006). 

log contains 3 x lO*" objects. The two catalogue utilize dif- 
ferent initial velocity distributions of NSs. One catalog (Ta- 
ble |2]l uses the initial velocity distribution of Arzoumanian et 
al. (2002), while the other (Table |3]l uses the initial velocity 
distribution of Faucher-Giguere & Kaspi (2006). Also derived 
are catalogue of simulated NS positions and proper motions 
with respect to an observer found at the solar circle (tables |4] 
and©. 

The space distribution at the cuiTent epoch obtained by the 
two initial velocity distributions implemented here, are some- 
what different. For example, I find that the resulting sky sur- 
face density, based on the two different initial velocity distri- 
butions, differs by up to 65%. The main differences between 
the two velocity distributions are in the low- and high-tails of 
the NSs velocity distribution (see Fig.lTJ. 

The space and velocity distributions of Galactic NSs were 
estimated by several past works. In Table |7] I compare some 
of the basic results of the simulations presented here (e.g., 
local density of NSs; scale height), with the ones obtained by 
previous efforts. 
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Fig. 10. — Sky surface density distribution of NSs, at tlie current epoch, for an observer located on the solar circle. The maps are presented using the Aitoff 
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